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ABSTRACT 

We have produced sensitive, high-resolution radio maps of 12 submillimetre (submm) galax- 
ies (SMGs) in the Lockman Hole using combined Multi-Element Radio-Linked Interferom- 
eter Network (MERLIN) and Very Large Array (VLA) data at a frequency of 1.4 GHz. Inte- 
grating for 350 hr yielded an r.m.s. noise of 6.0 /iJy beam^^ and a resolution of 0.2-0.5 arcsec. 
For the first time, wide-field data from the two arrays have been combined in the (u, v) plane 
and the bandwidth smearing response of the VLA data has been removed. All of the SMGs are 
detected in our maps as well as sources comprising a non-submm luminous control sample. 
We find evidence that SMGs are more extended than the general /iJy radio population and 
that therefore, unlike in local ultraluminous infrared galaxies (ULIRGs), the starburst compo- 
nent of the radio emission is extended and not confined to the galactic nucleus. For the eight 
sources with redshifts we measure linear sizes between 1 and 8kpc with a median of 5kpc. 
Therefore, they are in general larger than local ULIRGs which may support an early-stage 
merger scenario for the starburst trigger. X-rays betray active galactic nuclei (AGN) in six of 
the 33 sources in the combined sample. All but one of these are in the control sample, sug- 
gesting a lower incidence of AGN amongst the submm-luminous galaxies which is, in turn, 
consistent with increased X-ray absorption in these dust-obscured starbursts. Only one of our 
sources is resolved into multiple, distinct components with our high-resolution data. Finally, 
compared to a previous study of faint radio sources in the GOODS-N field we find systemati- 
cally smaller source sizes and no evidence for a tail extending to ~4 arcsec. Possible reasons 
for this are discussed. 
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1 INTRODUCTION 

SMGs form a class of object that was unknown before the ad- 
vent of modern bolometer arrays - starting with the Submillimetre 
Common User Bo lometer Arra y (SCUB A; Holland et al. 1999) - 
in the late 1990s jSmail. Ivison & Blainll 199711 . Effective follow- 
up of these sources has been particularly difficult due to the 
large positional uncertainty of the submm detections. As a result, 
most of our understanding of the properties of these sources has 
been made possible by deep, wide-field radio imaging, which ex- 
ploits the radio/far-infrared (far-IR) co rrelation to produce posi - 



nate predominantly from a starburst, there is likely to be a con- 
tribution from an AGN in some source s. Signatures of AGN 
include X-ray emission tAlexander et al.l 12003. 2005), compact 



ploits the radio/tar-mtrared (tar-lKj co rrelation to produce posi - 
tions with sub-arcsecond accuracy (e.g. I lvison et al.l 200ll2007h 



It is now k nown that SMGs lie at high redshift jChapman et al.l 
I2003ll2005h and it is believed that their large far-IR luminosities 
are due to the ultraviolet (UV) light produced by intense star for- 
mation (^IOOOMq yr~^) being reprocessed by dust; as such they 
are major contributors to the star-formation history of the Universe 
dHughes e t al. 1998). 

Although the emission from SMGs is believed to origi- 
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or je t /lobe radio emissio n I Chapman et alj 20041 ; IMuxIow et al.l 
l2005l ; lGarrettetal.ll200lh . unusual radio spectral ind ices as well 
as ce rtain line characteristics seen in optical |Swi nbank et id] 
'2004*) and mid-IR spectroscopy ^Menendez-Delmestre et al. 2007; 
iLutz et al.n2005i) . Star formation and AGN activity are believed 
to be closely linked due to the well-established correlation be- 
twe en black hole and bulge mass in nearby e l liptical galax- 
ies (lMagorrianetalj[T99i ; iMarconi & Hun3l2003l ; iHiiring & Rij3 
l2004h and, furthermore, feedback mechanisms have b een proposed 
whereby star formation is regulated by the AGN (e.g. Silk & Ree3 
1 19981 ; iDi Matteo, Springel & Hernquisll l2005t) . IPage et al.l j2004h 
have suggested an evolutionary sequence within which the tran- 
sition from X-ray absorbed to unabsorbed AGN represents the ter- 
mination of galaxy growth. 

In order to investigate the origin of the radio emission in 
SMGs, we have obtained high-resolution imaging of the Lock- 
man Hole with MERLIN, a grouping of seven telescopes across 
the United Kingdom. The maximum distance between telescopes, 
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Table 1. Lockman Hole SMGs included in our sample. Positions a nd 1.4-GHz flux densi ties are taken from the VLA catalogue of iBiggs & IvisonI j2006l 
Redshifts are determined spectroscopically and are taken from eiflier lChapman et alj j2005l) or llvison et all j2005ll2007l) . 



ID 


SHADES ID 


MAMBO ID 


Bolocam ID 


Radio position (J2000) 


Si A (^Jy) 


z 


SMGOl 


LOCK850.01 


LH 1200.005 


1 100.014 


10:52:01.250 +57:24:45.76 


73 


2.148 


SMG02 


LOCK850.04 


LH 1200.003 




10:52:04.226 +57:26:55.46 


66 


1.48 


SMG03 


LOCK850.09 


LH 1200.042 




10:52:15.634 +57:25:04.27 


58 


1.853 


SMG04 


LOCK850.16 


LH 1200.096 




10:51:51.690 +57:26:36.09 


110 


1.147 


SMG05 


LOCK850.17 


LH 1200.011 




10:51:58.018 +57:18:00.26 


93 


2.239 


SMG06 


LOCK850.30 






10:52:07.489 +57:19:04.01 


246 


2.689 


SMG07 


LOCK850.33 


LH 1200.012 




10:51:55.470 +57:23:12.75 


54 


2.686 


SMG08 


LOCK850.71 






10:52:19.086 +57:18:57.87 


98 




SMG09 


LOCK850.87 






10:51:53.362 +57:17:30.06 


93 




SMGIO 




LH 1200.008 




10:51:41.431 +57:19:51.90 


295 


1.212 


SMGU 






1100.003a 


10:52:13.376 +57:16:05.41 


249 




SMG12 






1100.0036 


10:52:12.256+57:15:49.59 


85 





215 km, corresponds to an angular resolution of ~200 milliarcsec 
(mas) at a frequency of 1.4 GHz, much higher than that ob- 
tained with the VLA at the same frequency. The Lockman 
Hole is notable for its very low Galactic Hi column density 
( iLockman. Jahoda & McCammonlll986) and has been targeted by 
a number of different instruments at diffe r ent wavelengths, in- 
cluding X MM-Newton /Hasinger et al."200l'; 'Mainieri et al."2002"; 
iBrunner et al. 2007), 5p/rzer ( Huang et al. 2004; Egami et al. 2004) 
and th e VLA dPe Ruiter et alJl997l : ICihegi et al.l2000l : llvison etal] 
l2002h . 

In this paper we combine our new MERLIN data with existing 
VLA data in order to increase the sensitivity to extended emission 
that would otherwise be resolved out or poorly constrained due to 
the lack of short spacings, whilst still retaining significantly higher 
resolution than can be obtained with the VLA alone. The paper is 
structured as follows. In the next section we describe the sample of 
SMGs as well as a control sample with which the structures of the 
SMGs can be compared. In f|3]we describe the observations and 
the data analysis, including a new technique for combining radio 
interferometric wide-field data in the Fourier plane, and present the 
maps. Discussion of our results and conclusions follow in 21™d[5] 



2 SAMPLE SELECTION 
2.1 SMGs 

In the submm waveband, the Lockman Hole has been targeted 
in several overlap ping surveys: at 850 fim jScott et ^ 1 20021 : 
ICoppin et all |2006|) with SCUBA on the James Clerk Maxwell 
Telescope (JCMT), at 1200/im with the Max-Planck Millimeter 
Bolometer array (MAMBO) on the Institu t de Radio Astron omic 
Millimetrique's (IRAM's) 30-m tel escope iOreve et a ljl2004l) and 
at 35 and 1100 p.m with SHARC II ^Laurent et al.l20 06) and Bolo- 
cam ( iLaurent et al.ll200^ , respectively, at the Caltech Submillime- 
ter Observatory (CSO). The Lockman Hole port ion of the SCUBA 
Half-Degree Extra galactic Survey (SHADES : iMortier et al.ll2005l : 
ICoppin et alj|200^ found the largest number of sources: 60 SMGs 
at 850 /xm over an area of approxim ately 360arcmin' ^ . This in- 
cluded data from the 8-mJy Survey ( IScottetalJl2002l :l Fox et al.l 
[2002). Combining source lists from several surveys h as been used 
as a check on the robustness of individual sources ilvison et al.l 
I2OO5I) . 



Our sample was compiled from the above references and com- 
prises all SMGs within 6 arcmin of the MERLIN and VLA pointing 
centre (lO'' 52™ 08':82, +57° 21' 33'.'8, J2000) that have a peak ra- 
dio flux density (as measured with the VLA) in excess of 50 ^Jy. 
The radial cutoff is required to ensure that the radio data are mini- 
mally effected by bandwidth and time smearing (i ]3.3t and primary 
beam effects. The primary beam for unequal apertures has been 
investigated by Strom (2004) and for baselines with the most dis- 
parate telescope sizes (76 m-25 m) is equal to 15.6 arcmin (FWHM), 
significantly in excess of our 12-arcmin diameter field. 

With these constraints we obtain a sample of 12 SMGs which 
we list in Table [Tj nine of which were detected by SHADES. A 
number of these sources are also detected by MAMBO and or 
Bolocam/SHARC II. The three remaining sources were only de- 
tected by either MAMBO or Bolocam/SHARC II; the radio iden- 
tification (ID) is unambiguous in e ach case. For LH 1200.008 
( iGreve et alj|2004l : llvison et alj |2005'), there is only a single, rel- 
atively bright radio source located near the su bmm position. The 
remaining SMG, 1 100.003 jLaurent et alj200^ . was originally de- 
tec ted by Bolocam and shown by later observations with SHARC 
II jLaurent et al.l 1200^ to be a blend of two sources. Both of 
these have a radio counterpart; we refer to them as 1 100.003a and 
1100.003b respectively, where 1100.003a is much the brighter of 
the two in the radio. Flux den s ities a nd positions for all sources 
are taken from lBiggs & IvisonI j2006l) a nd spectroscopic redshifts 
(available for seven of th e SMGs) from lchapman et all ( |2005|) or 
llvison et al. U2OO5II2OO7I) . 



2.2 Control sample 

As well as the SMG sample, we have defined a control sample with 
which the properties of the SMGs can be contrasted. This sample 
is made up of the remaining VLA detections within the 6-arcmin 
radius area that have peak brightnesses greater than 50 fiJy beam~^ 
i.e. it is defined in an identical fashion to the SMG sample. Of these, 
we exclude one obviously resolved source which is much brighter 
than the others, its total flux density exceeding 1 mjy. This gives a 
control sample of 21 sources. 
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Figure 1. VLA 1 .4-GHz map of the Lockman Hole, showing the positions of the SMGs (circles) and the control sample (squares). Also marked is the 6-arcmin 
radial cutoff used in defining our samples and the centre of the field. 



3 OBSERVATIONS AND DATA REDUCTION 



3.1 VLA 



The VLA data were first discussed in llvisonet aP ( l2002l) and we 
give only a brief summary here. The standard VLA wide-field 
'pseudo continuum' mode was used, two separate frequency bands 
near 1400 MHz (both containing right and left circular polarisa- 
tions, RCP and LCP) split into seven 3.125-MHz-wide channels 
each. The correlator integration time was 5 s. The full primary 
beam was mapped, with IMAGR, and the data self-calibrated us- 
ing a model of the sky (including bright sources in the primary 
beam sidelobes). The data from the several days of observations 
were then combined with STUFFR and produced final maps with an 
r.m.s. noise of approximately 5 /ijybeam^^. In order to simplify 
later mapping of the combined MERLINA^LA dataset, the VLA 
data had all CLEAN components which lay further than 6 arcmin 
from the phase centre subtracted using UVSUB. In this way, map- 
ping with the very large combined dataset could be confined to the 



central regions of the primary beam, the sidelobe response from 
brighter, more distant sources already having been removed. 



3.2 MERLIN 

The MERLIN data were observed over the period 2005 March to 
June using all seven telescopes, including the 76-m Lovell tele- 
scope which approximately doubles the sensitivity. A single 16- 
MHz band with a central frequency of 1408 MHz was split by 
the correlator into 32 separate 0.5-MHz wide channels with data 
produced every 4 s. The receivers were sensitive to both RCP and 
LCP, but a lack of correlator capacity does not allow the cross- 
correlation of the RCP and LCP data necessary for polarisation 
imaging. Phase referencing was accomplished using the nearby 
source, J1058-I-5628, with approximately 6 min spent observing the 
target and 2 min on the calibrator. In total, the Lockman Hole was 
observed for ~290h; this corresponds to an expected image sen- 
sitivity of ~6/iJy beam~^. Absolute flux calibration was based 
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Figure 3 - continued 



on observations of 3C286 and bandpass calibration with either 
OQ208 or DA 193. 



Initial flux and bandpass calibration took place in the usual 
manner with subsequent processing undertaken using the National 
Radio Astronomy Observatory's AIPS package. Gain solutions, first 
in phase and then in amplitude and phase were derived from the 
phase calibrator data and interpolated onto the target. The cali- 
brated data were subsequently combine d into a single dataset using 
STUFFR. As with lMuxlow et alj ( l2005h , four of the outer channels 
were not included when mapping as these reduced the image sen- 
sitivity. At this stage it is usual to scale the data weights to reflect 
the nominal sensitivity of each telescope. In preference to this, we 
have used FIXWT to set the weights based on the scatter in the data; 
this causes the weights to reflect the relative sensitivity of each tele- 



scope, additionally as a function of time (caused by, for example, 
the changing antenna elevation). 

In preparation for combining these data with those from the 
VLA, maps were made of sources identified with the VLA us- 
ing IMAGR with natural weighting. Those sources which were de- 
tected and were located further than 6 arcmin from the phase centre 
had their CLEAN components subtracted from the data, as with the 
VLA. Due to the different primary beam response of baselines that 
include the Lovell telescope, the MERLIN data were split in two, 
separating non-Lovell baselines from those that include the Lovell. 
Each subarray was then mapped separately and the distant sources 
subtracted from the data. A complete MERLIN dataset was then 
reconstituted by combining the two datasets with DBCON. 

Self-calibration of the data was attempted prior to component 
subtraction. Unfortunately, due to the lack of sufficiently bright 
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Table 2. Radio sources with peak flux densities greater than 50 fijy beam" ^ 
which lie within 6 arcmin of the phase centre and which do not have submm 
detections. These constitute our contr ol sample. Flux densitie s and positions 
are taken from the VLA catalogue o f iBiggs & Ivisoi] i200el) . 



ID 


Radio position (J2000) 


51.4 (^Jy) 


COl 


10:51:25.892 +57:21:54.61 


7S 


C02 


10:51:39.588 +57:17:32.06 


5S 


C03 


10:51:44.431 +57:22:03.74 


54 


C04 


10:51:44.492 +57:17:19.66 


137 


C05 


10:51:49.939 +57:25:21.73 


83 


C06 


10:51:50.113 +57:26:35.73 


115 


C07 


10:51:58.150+57:26:01.13 


62 


COS 


10:51:58.919 +57:23:30.13 


224 


C09 


10:52:05.491 +57:17:09.99 


85 


CIO 


10:52:12.493 +57:24:52.91 


222 


CU 


10:52:13.283 +57:26:50.36 


52 


C12 


10:52:16.903 +57:20:17.24 


lOS 


C13 


10:52:17.615 +57:21:27.37 


160 


C14 


10:52:29.740+57:17:53.31 


89 


C15 


10:52:31.522 +57:17:51.67 


56 


C16 


10:52:31.849 +57:24:50.39 


55 


C17 


10:52:32.877 +57:25:42.74 


115 


CIS 


10:52:33.748 +57:17:58.85 


62 


C19 


10:52:37.370+57:21:48.15 


155 


C20 


10:52:39.600+57:24:31.37 


143 


C21 


10:52:42.418 +57:19:14.59 


253 



sources, it proved impossible to find enough solutions with suffi- 
cient signal-to-noise over the short timescales that were hkely to 
characterise the residual phase errors. However, the phase and am- 
phtude solutions are generally smooth functions of time (the obser- 
vations having been made in good weather and close to solar min- 
imum) and the residual phase variations are likely to be small. At 
lower elevations the relative phase varies more rapidly due to the 
greater differential path length, but these data have lower weight 
and thus contribute less to the images. We also note that the angle 
between the calibrator and target field is small, 1?2. 

3.3 Bandwidth and time smearing 

A problem inherent in imaging wide fields with radio interferom- 
eters is chromatic aberration caused by the finite width of the fre- 
quency channels. This results in a source being stretched along the 
vector pointing from the source to the phase centre (i.e. radially) 
with a magnitude proportional to a factor equal to the distance from 
the phase centre multiplied by the fractional bandwidth. For our 
data, the effect has been minimised by splitting the bandwidth into 
narrow channels and by restricting our sample to sources within 
6 arcmin of the phase centre. 

The apparent size of a smeared source is a convolution of the 
radial smearing factor with the telescope beam; therefore smeared 
sources appear to be much bigger when observed with the VLA 
than with MERLIN, even if the smearing factor is the same. In 
order to remove the remaining VLA smearing we have replaced the 
VLA data corresponding to each source with a de-smeared version. 
This is possible because at the resolution of the VLA all our sources 
have very simple structures and can be accurately modelled with 
single Gaussian components. These can be subtracted from the data 
and replaced with new model components that have been corrected 
for the appropriate amount of smearing. 



The Gaussian model fitting is performed on the VLA maps us- 
ing the task JMFIT. This has a parameter (BWSMEAR) which, when 
set to the fractional bandwidth of the observations, corrects the fits 
for bandwidth smearing and deconvolves these with the smeared 
beam. This gives the source size that the VLA would have mea- 
sured in the absence of smearing. With the fit parameters fixed at 
the values found in the first run, JMFIT is run again, this time with- 
out the smearing correction. This causes the deconvolution to be 
performed with the (unsmeared) CLEAN beam and thus returns the 
size of the source as contained in the data i.e. the smeared size. 
The smeared and unsmeared model components are then subtracted 
from and added to the data respectively, using UVSUB. 

Averaging visibility data in the time domain causes a simi- 
lar distortion to bandwidth smearing, but in the opposite direction 
i.e. tangentially. Sources have visibility functions which vary more 
rapidly with time at increasing distance from the phase centre and 
so averaging restrictions are accordingly tighter. With similar cor- 
relator averaging times, the MERLIN data will have the more re- 
stricted field of view of the two arrays. However, for our 6-arcmin 
field of view we find that the effect of time smearing is negligible 
for the VLA and w ill reduce the visibility of a MERLIN baseline 
by only 1 per cent jThompson. Moran & Swensonlll98^ . The ef- 
fect in the combined maps is also negligible given that the restoring 
beam is twice that of MERLIN alone. 

3.4 Combining the MERLIN and VLA data 

Combining data from different sized arrays is a reasonably stan- 
dard procedure in radio astronomy, the resulting dataset contain- 
ing a greater range of projected baseline spacings and hence giving 
simultaneous sensitivity to a wider range of source sizes. This is 
perhaps most commonly associated with the different array config- 
urations of the VLA, but the Australian Compact Telescope Array 
(ATCA) and the Westerbork Synthesis Radio Telescope (WSRT) 
also routinely combine data from separate configurations. The tech- 
nique has also been used to combine data from completely dif- 
ferent instruments, to produce hybrid MERLIN/VLA and MER- 
LIN/European VLBI Network (EVN) datasets, for example. 

Whilst the combination can be accomplished in the image 
plane, it is instead usually done in the Fourier plane. This is be- 
cause the improved (u, v) coverage is better able to constrain de- 
convolution algorithms such as CLEAN and which should there- 
fore produce more reliable images. Having access to the visibili- 
ties has a further advantage in that the Cotton-Schwab CLEAN al- 
gorithm (Schwab 1984) can be used (as implemented in IMAGR). 
This is superior to image-plane deconvolution (such as is offered 
by APCLN in AIPS) for a variety of reasons, including more accu- 
rate subtraction of the CLEAN components and the ability to map 
multiple fields simultaneously. 

Whilst relatively straightforward with the single-channel 
datasets that are used in the majority of radio continuum mapping, 
the situation is more complicated when one is trying to image wide- 
field data containing multiple channels. Within AIPS, the multiple 
channels are grouped together in so-called IFs (the AIPS name for a 
frequency band) with the (u, v) coordinate of a visibility in a given 
channel calculable from the reference frequency of the IF and the 
channel width, information which is available from the header of 
the dataset. As these numbers are generally different for a VLA and 
MERLIN multi-channel dataset, it is impossible to simply combine 
each by merging the IFs (as is usually done using DBCON); calcu- 
lating the correct («, v) cooordinate o f a visibility would b e im- 
possible. The problem is discussed by IMuxIow et ai] ( l2005h who 
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Figure 4. Difference in position angle of the source radial direction (from 
the pliase centre) and tlie measured source major axis as a function of dis- 
tance from the phase centre. 



instead chose to combine the Dirty Images and deconvolve using 
APCLN. 

It is, however, possible to combine the data in the Fourier 
plane, by converting a multi-channel dataset into one with only a 
single channel. This is achieved by extracting each channel from 
its IF using tJVCOP, and then re-combining with DBCON. During 
this last stage the (u, v) coordinates of each channel's data are set 
to their correct, absolute, values and the resultant dataset has the 
same structure as one with a single channel. After performing this 
separately on the two IFs of the VLA and the single IF of MER- 
LIN, the three datasets can be further combined into one. The tech- 
nique is conceptually simple and easily implemented in MPS. No 
re-weighting of data was done in DBCON as the weights of the 
MERLIN and VLA datasets were very similar. 



3.5 Mapping the combined data 

The combined MERLIN and VLA data were mapped using IMAGR. 
We used a robust weighting scheme, combining ROBtJST = with a 
Gaussian taper in the (u, v) plane. This produced a roughly circular 
beam (521 x 481 mas^ at a position angle of 23? 6); this corresponds 
to an improvement in resolution by a factor of nearly three relative 
to the VLA-only maps (or an order of magnitude in terms of beam 
area). Each map comprised 1024 x 1024 50-mas pixels and after 
deconvolution had an r.m.s. noise of ~6.0/iJy beam~^. We note 
that the quoted sensit ivity of the combined MERLINA'LA maps of 
lMuxloweta"l]j2005l) was better by a factor of two compared to the 
single-array maps. Despite this significant improvement, no new 
sources were found within the central 3-arcmin field. 

We also tried making maps with a smaller beam, approxi- 
mately equal to the MERLIN-only beam (~200 mas). However, in 
order to achieve such a beamsize, the VLA data must necessarily 
be down-weighted and thus contribute less to the map. These maps 
looked very similar to the MERLIN-only maps and where there is 
extended emission it is generally weak, patchy and over-resolved. 
We note that the ma jority of the sources (61 per cent) mapped by 
iMuxlow et 'Z H2005t) were restored with a 500-mas circular beam, 
similar to our own. 



MERLIN/ VLA maps with an angular extent of 7.25 x 
7.25 arcsec^ of the twelve SMGs in our sample are shown in Fig.[2] 
Contours are plotted at —1, I, 2, 4, 8, etc. times the 3 a noise and 
greyscales are shown in units of /ijybeam^^. Maps of sources in 
the control sample can be seen in Fig. [3] Clear detections are made 
of all the SMGs and control sources. 



4 DISCUSSION 

The sources in our combined sample (SMGs and control sources) 
display a variety of morphologies, but we first note that our ~0.5- 
arcsec beam has not resolved any of the sources into separate, mul- 
tiple components although most are extended. In order to quantify 
some properties of the sources, we have fitted a single Gaussian to 
each using JMFIT. Where the sources are weak or show evidence 
of multiple components the fits will be less reliable, but will still 
allow us to measure the approximate source size and orientation. 

Imaging such large fields with radio interferometers can lead 
to the suspicion that the observed source structures are being af- 
fected by smearing. We are confident, however, that the steps we 
have taken (restricting the sources that we image to those within 
a 6-arcmin radius and correcting the VLA data for the effects of 
bandwidth smearing) have resulted in our maps being unaffected 
by this effect to any great degree. To test this we have computed the 
difference between the measured position angle of the major axis 
of the Gaussian fit and the direction of the source from the field 
centre. Fig. |4] shows the difference in these angles as a function of 
radial distance; the signature of smearing, i.e. increased clustering 
around zero or 90 degrees towards larger radial distances, is not 
seen. 



4.1 Size of the radio-emitting region of SMGs 

The deconvolved sizes (FWHM) of the fitted Gaussians are shown in 
Table[3]along with a measure of the relative goodness of fit (the sum 
of the square of the residuals between the model and the data). The 
distribution of this un-normalised chi-square statistic has an aver- 
age of 3. 7 ±0.8 (x 10~* Jy^ beam^"^) with two significant outliers, 
SMG 1 and C 1 3 . An examination of the maps for these two sources 
shows that each is significantly resolved; the fitted sizes for these 
sources should be viewed with caution. Also included in Table [3] 
are the errors on the deconvolved sizes determined by JMFIT. 

Based on the JMFIT error analysis, all of the sources are re- 
solved, with deconvolved sizes greater than the associated I-cr er- 
ror. Three sources, however, have a major axis and uncertainty that 
are of similar magnitude (within a factor of two) and a minimum 
size for their major axes (as returned by JMFIT) of zero (source 
smaller than the beam). These three sources are the least likely to 
have been resolved and they (C06, C 1 1 and C 1 8) are all from the 
control sample. 

In Fig. |5] we plot the deconvolved major axes for our sam- 
ple of SMGs. These cover the range 0.2-1.0 arcsec, with a slight 
tendency for the larger sources to be more common. The me- 
dian diameter is equal to ■~0.65±0.1 arcsec. A previous investi- 
ga tion into the size of th e radio emission from SMGs is reported 
bv lChapman et alj ( l2004l) . using the combined MERLINA'LA data 
from GOODS-N that is presented in more detail bv lMuxlow et al.l 
( I2OO5I) . [chapman et al. measure a median source diameter of 
0.83 ± 0.14 arcsec, larger although not incompatible with our own 
measurement. 

A comparison of the two values is complicated by the different 
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Table 3. Deconvolved source sizes for each source. The relative goodness 
of fit is indicated by the sum of squares of residuals. 



Source 


Major axis 
(arcsec) 


Minor axis 
(arcsec) 


Sum of square of residuals 
(xlO-* Jy^beam-^) 


SMGOl 


0.55 ± 0.10 


0.29 ± 0.14 


3.4 


SMG02 


0.84 ± 0.14 


0.31 ± 0.16 


3.6 


SMG03 


0.98 ± 0.32 


0.81 ± 0.31 


3.0 


SMG04 


0.92 ± 0.12 


0.59 ± 0.10 


2.9 


SMG()5 


0.94 ± 0.14 


0.71 ± 0.13 


3.6 


SMG()6 


0.19 ± 0.05 


0.13 ± 0.08 


2.9 


SMG07 


0.39 ± 0.14 


0.20 ± 0.19 


4.3 


SMG()8 


0.40 ± 0.09 


0.30 ± 0.10 


4.4 


SMG09 


0.80 ± 0.13 


0.41 ± 0.13 


4.2 


SMGIO 


0.61 ± 0.04 


0.44 ± 0.03 


9.7 


SMGll 


0.34 ± 0.04 


0.22 ± 0.06 


2.0 


SMG12 


0.66 ± 0.11 


0.40 ± 0.12 


3.0 


COl 


1.16 ± 0.20 


0.40 ± 0.17 


2.9 


C02 


0.52 ± 0.11 


0.10 ± 0.14 


2.9 


C03 


0.49 ± 0.13 


0.41 ± 0.14 


5.4 


C04 


0.36 ± 0.06 


0.16 ± 0.11 


4.3 


C05 


0.46 ± 0.09 


0.33 ± 0.09 


3.6 


C06 


0.14 ± 0.12 


0.03 ± 0.10 


3.3 


C07 


0.83 ± 0.16 


0.58 ± 0.15 


3.9 


COS 


0.20 ± 0.05 


0.00 ± 0.06 


4.8 


C09 


0.43 ± 0.15 


0.36 ± 0.16 


3.9 


COlO 


0.30 ± 0.06 


0.25 ± 0.07 


3.9 


coil 


0.25 ± 0.19 


0.12 ± 0.15 


4.0 


C012 


0.36 ± 0.08 


0.25 ± 0.10 


3.1 


C013 


0.66 ± 0.08 


0.55 ±0.08 


8.0 


COW 


0.53 ± 0.09 


0.33 ±0.09 


4.4 


C015 


0.50 ±0.11 


0.05 ±0.12 


2.2 


C016 


1.39 ± 0.28 


0.38 ±0.24 


4.2 


C017 


0.52 ± 0.07 


0.09 ±0.11 


2.5 


COl 8 


0.26 ± 0.19 


0.00 ±0.10 


2.6 


C019 


1.29 ± 0.14 


0.37 ±0.10 


5.0 


C020 


0.36 ± 0.07 


0.19 ±0.14 


4.5 


C021 


1.09 ± 0.08 


0.76 ±0.07 


4.6 



ways in which the source sizes have been measured and the diffi- 
culty of measuring the sizes of weak, moderately resolved sources. 
Whilst we have quoted deconvolved major axes of a Gaussian fit, 
IChapman et al.l give the largest extent bounded by the 3ct contour. 
Gaussian fitting allows us to deconvolve the beam from the fit- 
ted component, something that is particularly important when the 
source sizes are of order the beam or smaller. Weak extensions are 
not well modelled, but contribute very little to the total flux density 
and would give unreliable size estimates in any case because of the 
low signal to noise. 

Measuring the sizes of SMGs is crucial to understanding their 
formation. In the local universe, ULI RGs seem to be analogue s 
of the high-redshift SMG population jSanders & Mirabe]||T99^ . 
The most luminous of these are predominantly merger events 
and there is evidence that the similarly luminous high-redshift 
SMGs are also commonly undergoing mergers; physical associa- 
tions are directly inferred from spectroscopy and multi-wavelength 



imagmg jlvison et alJI 19981: ISmail et alJI 19981: llvison et alj 



IConselice et"al] l2003l : ISmail et al.l |2004| : ISwinbank etaH 



20001 : 



2004h 



whilst the double-peaked profiles often seen i n mm-wave spec- 
troscopy of the vast reservoirs of molecul ar gas jFraver etal.l 19981 
ll999l : lNeri et al.ll2003l : ICreve et al.ll2005h are suggestive of orbital 
motions. 




Source size (arcsec) 

Figure 5. Histogram of the source sizes (deconvolved major axis) of the 
SMGs (solid) and control sample (cross-hatched) in the Lockman Hole. 



iTacconi et alj ( |2006[) present observations at 1 and 3 mm of 
SMGs with the IRAM Plateau de Bure interferometer which at the 
shortest wavelength have a beamsize (~0.6 arcsec) that is similar 
to that of our MERLIN/VLA maps. They find that the SMGs (as 
seen in CO and/or mm-wave continuum emission) are generally 
compact; their sample of ten sourc es has a media n size ^0.5 arcsec 
(±0.1 arcsec). The sizes quoted bv lTacconi et al.l are 'angular aver- 
aged' and to facilitate a proper comparison with our data we plot 
the average of the major and minor axes of our SMGs in Fig. |6] 
along with the' Tacconi et al .' results. The two distributions appear 
similar, in both the position of the peaks and the range of sizes cov- 
ered. A Kolmogorov-Smimov (K-S) test ( Press et al.lll992) gives 
a probability of 84 per cent that the differences between the two 
could arise by chance and thus strongly favours the null hypothesis 
that the two distributions are the same. This mirrors what is seen at 
much lower redshif t, where the radio emission (Condon et al. 199^ 
and molecular gas jPownes & SolomorJI 19981) have similar sizes. 

Converting to linear distances for the SMGs with redshift 
we fin d a range 1 .2-7.7 kpc with a median of 4.9 kDc: lTacconi et al. 
12006) measure a smaller, but similar, linear extent of <4kpc. 
Given that the size of the central starburst in local ULIRGs is 
< 1 kpc, the high-redshift starbursts appear to be in general signifi- 
cantly larger, although the smallest (SMG06) is of comparable size 
to th ose seen at low red shift. 

llvison et alj ( |2007|) have explored the fraction of SMGs that 
contain multiple radio counterparts and find evidence that this is 
more common for the brighter sources, with separations between 2 
and 6 arcsec0. They are unable to probe below 2 arcsec due to the 
resolution of the VLA-only data, but despite the smaller beam of 
our combined MERLIN A'LA maps we see no examples of multiple 
sources with separations <2 arcsec. We can increase the resolution 
still further, at the expense of accurate imaging of extended emis- 



0.27, n 



A 



^ We assume a flat 'concordance' cosmology with Hm 
0.73 and Ho = 71 kms^i Mpc^^ (Spergel et al. 2003). 
^ SMG02 (Lock 850.04) has two radio components close to the submm po- 
sition and although only one of these is included in our sample, the weaker 
of the two (<50^Jy) is detected in the MERLINA'LA map and is visible as 
the 3(T contour north-northeast of the stronger detection. 
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Source size (orcsec) 



Figure 6. A comparison of the sizes of SMGs as measured by our MER- 
LINA'LA radio data (solid) and bv .Tacconi et al., (.2006,) in the mm wave- 
band (cross-hatched). 



50 100 




3 2 10-1-2-3 

ARC SEC 



Figure 7. MERLIN-only map of SMG09 (LH 1200.008). The beam is 
shown in the bottom-left comer and has dimensions of 373 X 223 mas^ at a 
position angle of 9? 7. The greyscale has units of /jjy beam~^ and contours 
are drawn at —1, 1, 2, 4, 8, etc. times the 3 a noise (3 X 8.5 /ijy beam~^). 



sion, by examining the MERLIN-only maps. One of these, SMGIO, 
reveals that a prominent extension in the MERLIN/VLA map con- 
tains a relatively high signal-to-noise (~6) component (Fig. |7}. A 
two-Gaussian fit in JMFIT gives a separation of 0.6 arcsec and a flux 
density ratio of ~5:1. At the redshift of 1.212 this angular separa- 
tion corresponds to a linear distance of ~5 kpc. The nature of this 
second component is unknown and may be a radio jet, but it may 
also be a small-separation merger. 



o 




0.2 0.4 o.e 0.8 

Peak flux density ratio (MERLIN/VLA) 

Figure 8. Histogram of the peak flux density ratio (MERLIN/VLA) of the 
SMGs (solid) and control sample (cross-hatched) in the Lockman Hole. 

4.2 SMG sizes relative to the /ijy radio population 

Also plotted in Fig. |5] are the major axes of the fits to the con- 
trol sample, i.e. those ^ly radio sources which are not luminous in 
the submm waveband. The range of angular sizes for both samples 
is similar, covering approximately 0.1-1.4 arcsec although the me- 
dian size of the control sample is formally smaller than that of the 
SMGs, 0.5 arcsec compared to 0.65 arcsec (with 1 o errors of about 
0.1 arcsec). Another way of assessing the extension or otherwise of 
the two source populations is to compare the peak flux densities 
in the MERLIN-only and combined MERLIN/VLA maps. These 
ratios are plotted in Fig. [8] and a similar picture emerges; the con- 
trol source distribution contains sources which are, in general, more 
compact. A K-S analysis gives low probabilities that the distribu- 
tions are the same (30 and 17 per cent for the size and peak flux 
density plots, respectively). 

We plot the sizes of all our sources in Fig. |9] and include 
the bright source that was excluded from our control sample (a 
large, ~4 arcsec, source which contains a radio core and two-sided 
jet). The distribution peaks strongly at angular exte nts in the range 
0.4-0. 6 arc sec. This is in ap parent disagreement with lMuxlow et aU 
( I2005I) and lChapman et aljr2004l) who find that the average source 
size of /iJy radio sources is ~1 a rcsec and that th ere is a prominent 
tail extending to -^4 arcsec. The IMuxIow et al.l data are also plot- 
ted in Fig. |9]so that the two samples can be easily compared; the 
two distributions are very different and they can formally be re- 
jected as representing the same population with a K-S probability 
in excess of 99 per cent. 

The distributions in Fig.|9]differ in two ways. Firstly, the Lock- 
man Hole data appear to consist of only a single source popula- 
tion, this being distributed approximately normally. The GOODS- 
N distribution is similar for the smaller sources, but peaks at a 
larger source size and has a significant tail out to angular extents of 
~4 arcsec. The absence of this tail in our sample (save for a single 
source) is particularly striking. Even taking into account possible 

^ Our histogram of the IMuxIow et al.l data is not identical to that in the 
original publication, probably due to slight differences in the way the data 
are binned. 



© 0000 RAS, MNRAS 000, 000-000 



High-resolution radio observations ofSMGs 1 1 



o 




2 3 4 5 

Source size (arcsec) 



Figure 9. Histogram of the source sizes of our Lockman Hole sample as 
well as the faint radio population in the GOODS-N field of Muxlow et aO 
I2OO5 ). For the Lockman Hole (solid area) the sizes refer to the deconvolved 
major axis of a Gaussian fit whilst for GOODS-N (cross-hatched) it is the 
'largest angular size'. 



errors in our model fitting (relatively weak extensions are not well 
modelled), we can be certain that only one source in our 12-arcmin- 
diameter field is larger than 2 arcsec. 

The explanation for the difference in source size distributions 
in Fig. 121 especially with regard to the tail, is prob ably a combi- 
nation of the larger sample and lower flux limit of IMuxIow et al.l 
(40 ^Jy compared to our 50 ^Jy) and cosmic variance. Our sources 
may also be smaller, particularly for sources in the population 
which peaks at around 0.5-1 arcsec, because of differences in the 
way the sizes were measured and also because we have attempted 
to remove the bandwidth smearing from our VLA data, something 
that was not done for the GOODS-N data. 



5 CONCLUSIONS 

We have presented sensitive maps (6/iJybeam~^ r.m.s.) of 12 
SMGs in the Lockman Hole made from combined MERLIN and 
VLA data. This is the first time that wide-field (multi-channel) data 
from MERLIN and the VLA have been combined in the (u, v) 
plane. The procedure to do this is straightforward, practical with 
modem desktop computing facilities and can be accomplished with 
existing software tools, i.e. AIPS. 

From the radio morphology alone it has generally not been 
possible to determine unambiguously whether a source contains an 
AGN, except in one case which has a prominent jet (C19). In the lo- 
cal Universe, radio emission from most starburst galaxies is known 
to be concentrated in their galactic nuclei, so source compactness 
is no guarantee of an AGN. Indeed, even higher resolution - very 
long baseline interferometry - is required to improve our diagnos- 
tic ability, via sensitivity to more compact jets and via the ability 
to constrain the brightness temperature of unresolved emission to a 
regime unique to AGN. 

Our maps have resolved the vast majority of the SMGs as well 
as a complementary sample of radio sources of similar flux den- 
sity. A comparison of the properties of these two populations shows 



that the size distributions overlap considerably, but the SMGs are in 
general larger than the control sample (six of which harbour AGN, 
being luminous X-ray sources - see below). Even if some of our 
submm-luminous sources also contain significant emission from 
AGN, this suggests that the starburst component in high-redshift 
SMGs is extended and not confined to the nucleus. 

This conclusion is strengthened by the measurement of simi- 
lar sizes for the radio-emitting regions of SMGs and their molec- 
ular gas reservoirs, suggesting that AGN are not significantly con- 
tributing to the measured size of the radio emission. However, most 
high-redshift SMGs have much larger sizes (1.2-7.7 kpc) than local 
ULIRGs (^ 1 kpc), possibly the result of witnessing them at a less 
advanced stage of the merger process. 

Despite our high resolution, as compared with the VLA, we 
find only one source that is resolved into multiple, distinct compo- 
nents. Determining the true nature of this source (and SMG06) - 
merger-induced starburst or AGN plus jet? - will require follow- 
up observations, but small separation (<2 arcsec) multiple counter- 
parts to SMGs appear to be rare. 

X-ray emission is a good indicator of AGN activity. Based 
on co unterparts in the XMM-Newton catalogue of iBrunner et al.l 
( |2007|) . we find a lower detection rate of X-ray luminous AGN 
amongst the submm-luminous sources than amongst the control 
sample (Poissonian probability, 89 per cent). Six of the control 
sources have a likely X-ray counterpart (C05, CIO, C12, C13, C20 
and C21) versus only one of the submm galaxies (SMG06). SMG06 
is the hardest of the seven X-ray detections, consistent with en- 
hanced X-ray absorption in a dust-obscured starburst ( Pa ge et^ 
1 2004) . It is also the most c ompact radio emitt er amongst the SMGs, 
with a flat radio spectrum jlvison et alj2007b . 

e-MERLIN and the EVLA will soon be available, the broad- 
banded and hence much more sensitive versions of the current ar- 
rays. This will enable the /ijy radio population to be imaged at high 
resolution over larger areas and at increasingly faint flux limits. 
Also to come is the eSMA, a collaboration between the JCMT, the 
CSO and the telescopes of the Submillimeter Array. The resulting 
instrument will have a resolution at 345 GHz of 300 mas, extremely 
complementary to the angular scales being probed by the combined 
MERLIN and VLA maps presented in this paper. This will allow 
a direct comparison between the submm and radio morphologies 
of SMGs and we can thus examine the radio/far-IR correlation at 
unprecedented resolution. 
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